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ABSTRACT 


The primary objective of this project is to estimate the effectiveness of decoys 
against incoming missile attacks in a multiple-ship battle group ( BG ) environment. 

In this problem, we allow the number of ships, number of decoys, and the values 
of ships to vary. We also vary the values of the hard kill probabilities, the splash prob- 
abilities of decoyed missiles, probabilities that a missile’s lock is broken by seduction 
decovs and the quality factor of the distraction decovs. 

When an ASCM attacks a BG, it may hit a ship, it may get shot down, or it may 
get diverted. If it gets diverted it may lock onto a neighboring friendly ship. The 
measure of effectiveness ( MOE ) is the probability that all ships survive the missile at- 
tack. 
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I. INTRODUCTION 


A. THE PROBLEM 

The problem considered here is to defend a multiple-ship Battle Group ( BG ) 
against anti-ship cruise missile ( ASCM ) attacks. For this reason, a multiple-ship mis- 
sile attack model has been developed to study the effect of Electronic Warfare ( EW ) 
expendable decoys, which may increase the survivability of a single ship or a battle 
group. The effectiveness of a single decoy, or of a group of identical or different decoys, 
employed during an ASCM attack on a BG may be degraded as a result of missile di- 
version onto neighboring own-force or allied ships in the BG. 

The model that we present in Chapter ITI is intended to investigate how expendable 
decoys may increase survivability. The emphasis is on simplicity consistent with repres- 
enting the multi-ship situation rather than on complete realism. Individual systems may 
have peculiarities that may make some of the assumptions in the model invalid. 

In this study it 1s assumed that the BG is being attacked by one or more guided 


missiles. For each individual missile, an attack is assumed to consist of two phases: 
1. Target acquisition by the missile seeker or guidance system. 


2. Missile homing to impact. 


For an individual missile to successfully engage a target, it must complete both 
phases. Expendable decoys can disrupt either of these phases. If the missile cannot be 
distracted during the target acquisition phase then it might be ” seduced ” in its homing 
process ( breaking the acquisition of the missile seeker ). The missile might then be di- 
verted onto another ship, or an empty zone ( occupied by no ships), or the diverted 
missile might ” splash-down ” to the ocean ( missile might divert not only from its course 
but also it might divert to an altitude closer to the ocean level ). Hence the missile might 
“ splash “ into the ocean as a result of losing its altitude after seduction. Also, the ship 


air defense guns might shoot down the missile. If the missile is not seduced then it might 


be shot down by ship air defense guns, or if it is not shot down by ship guns it may hit 
a ship in the BG. 

For many missiles vs many targets ( ships ), it is assumed that the missiles act in- 
dependently and their errors are not correlated. It is also assumed that the effects of the 
decoys are independent. These assumptions of independence may hold in some opera- 
tional circumstances but in others they may not hold. 

The methods of defense against ASCM attacks for a BG that has at least one or 
two high-value ships and various combat and escort ships taken into account are: 

1. Hard kill by ship guns or antimussile missiles. 


2. Soft kill by using EW assets ( splash down or diversion by using decoys ). 


In this study, we developed a simple model that enables us to compare different 
factors, such as effectiveness of different types of decoys ( seduction and distraction ), 
the effect of changing the numbers and types of decoys, decoys vs hard kill, and non- 
linear interactions caused by diversion onto friendly ships in a multiple-ship BG. 

In this study, the following conditions hold for ASCM attacks on a multiple-ship 
BG environment situation: 


1. Effectiveness is measured as survivability of the entire BG; not only for a single 
ship. 


2. Multiple-ship environment; one or two-high value ships such as aircraft carriers, 
battle ships, or supply ships in or near the center of the BG. Low-value ships, such 
as destroyers, are far from the center of the BG. 


3. Nonlinearities; ship-missile-ship interactions are important, not just one ship- 
mussile/mussiles interaction. 


4. Decoys can be expendable decoys deployed from the ships for defending the ship 
itself, or the BG, against incoming ASCM attacks. Two types of expendable decoys 
are considered: 


a. Distraction decovs provide alternate targets for the missile during the target 
acquisition process. 


b. Seduction decoys disrupt the homing process. 


Since this 1s an unclassified study, some additional assumptions regarding perform- 
ance parameters are made to exercise the model. The model, computer code, and sensi- 
tivity studies can be readjusted for special conditions. Hence, the model and sensitivity 
studies give results that depend on unclassified information. 


Il. BACKGROUND 


A. BATTLE GROUP (BG) 
As depicted in Figure 1, an idealized BG is a group of ships that has at least one 


or two high-value ships and other low-value ships. 





Figure 1. An idealized battle proup 


The approximate number of total ships in the group is 8-15, which can be divided into 


two categories of ships based on their importance. These categories are: 
1. High-value ships 


2. Low-value ships 


High value ships may be aircraft carriers (CV's), Battle ships (BB’s), and supply 
ships (AOE’s). Low value ships may be guided missile cruisers (CG’s), destroyers (DD's), 
and frigates (FF’s). The values of the ships in a BG may be related to the mission of the 
BG , location of the operation scene, and other factors that affect the BG scenario. 

The basic weapons of a BG for defending itself against ASCM attacks are: 

I. Guns. 
2. Antimissile missiles. 


3. Electronic Warfare assets, such as electronic counter measures 


( ECM ), electronic counter counter measures (ECCM), and decoys for seduction or 
distraction. 

In this study we are primarily interested in the weapons for hard kills ( guns, anti- 
missile missiles ) and decoys for soft kills. 

Mission types for a BG usually require it to remain on the ocean for a long period 
of time and to conduct some independent operations. Generally, the main threats to a 
BG come from the air, under the water, and the surface. But, in modern warfare, the 
most likely unexpected method of attacking a BG can be from under water or air. Today 
most conventional submarines may use ASCM’s to attack surface targets. Hence we 


discuss ASCM attacks from different platforms. 


B. ANYISHIP CRUISE MISSILES ( ASCM’S ) 

ASCM attacks on a BG may come from any direction relative to the course of the 
BG. ASCM’‘s can follow different flight patterns. Two possible patterns can be seen in 
Figure 2. The processes that an ASCM goes through in its flight pattern are: 

1. Target acquisition process. 


2. Homing process to impact. 


The target acquisition process for an ASCM includes tracking, locating, identifica- 
tion and classification of the target; calculations for fire control; launching the missile; 
flight until the final phase, which includes precise tracking to the prelocated position of 
the target ( homing to impact ). The target acquisition process can vary for different 
tvpes of missiles. We assume that this process occurs approximately from launching to 
6000 - 20000 yards from the target that the anti-ship cruise missile ( ASCM ) 1s aimed 
at. The range of this process for a missile depends upon the characteristics of the guid- 


ance svstem of the missile. 
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Figure 2. Target acquisition and homing processes of different ASCM’s 


The homing process to impact for an ASCM starts at the end of the target acquisi- 
tion process and continues to impact or muss of the target. The homung process occurs 
in relatively short distances compared to the target acquisition process. Although it 
varies from missile to mussile, we assume that homing process can occur around 6000 - 
20000 yards. 

The speed of the incoming missile can vary from mach 1.5 to mach 2.5 for different 


tvpes of missiles. 


C. INTERACTIONS, ASSUMPTIONS, AND THE NATURE OF THE PROBLEM 
The speed of the BG is usually from 10 to 30 knots. Considering the area that is 
covered by a BG, and comparing the speeds of the BG and the ASCM, we can assume 


that because the speed factor of the BG is small compared to the speed factor of the 


incoming mussiles, speed does not have a major effect on avoidance maneuvers of the 
BG against ASCM’s. For this reason we can assume that the BG is almost fixed during 
the attack. Hence, we consider that the distance traveled by the ASCM in the target 
acquisition process is longer than the distance traveled by the ASCM in the homing 
process. Thus, this assumption about the fixed location of the BG is even more justified 
in the homing process than the target acquisition process. 

As previously mentioned, a BG that has different kinds and values of ships can be 
attacked by ASCM’s launched by surface , air or underwater platforms. In each case 
ASCM’s are the faster, most unexpected, precise, and preferred manner of attack against 
a BG. Each different type of missile usually shows almost the same flight characteristics. 

Hard kill devices of the BG for defending itself against ASCM’s are high speed small 
caliber air defense guns, such as 20 millimeters, 35 millimeters and 40 millimeters --- for 
example, the PHALANX system. Short range surface to air missiles ( SAM’s ) are also 
used --- for example the SEA SPARROW air defense missile system. Since this study 
contains unclassified information, a hypothetical hard kill probability ( PHK ) of the BG 
for an incoming missile is chosen. ; 

Soft kill devices ( decoys ) are of two types: distraction decovs and seduction decoys. 
The quality factor ( Q ) for distraction decoys can be defined as a measurement of the 
distraction abilitv of a distraction decoy. If a distraction decoy has Q = 1.0, it has the 
same effect as a real ship on a missile e.g. essentially the same Radar Cross Section ( 
RCS ). Quality factors take on values in the interval [0,1] 

The probability of breaking the lock of a missile by a seduction decoy, ( B ), can 
be defined as a measurement of the seduction ability of a seduction decoy. After se- 
ducing an incoming missile, if the seduction decoy successfully breaks the lock of the 
incoming missile, the missile may divert or splash into the ocean. For different flight 
patterns there are different splash probabilities. For example, the flight pattern in Fig. 
2.a has a higher splash probability than the one in Fig. 2.b 


Il. SURVIVABILITY OF BATTLE GROUP 


A. THE MODEL TO CALCULATE THE PROBABILITY MASS FUNCTION 
VALUES ( PMF ) FOR NUMBER OF ENCOUNTERS. 


An idealized BG, hes approximately 8 - 15 ships. If we define N,, A, and, R by 
1. N,: number of ships in the BG, 
2. A: area covered by the BG, 
3. R: radius of the BG, 


then the ship density (p) for the BG which covers area A 1s 


It seems reasonable to suppose the ships will maintain their distance from each other 
when an attack is imminent. Hence for simplicity of analysis we assume the ships are 
located on the intersections of the grid lines within the area A. Later we will allow the 
ships to be uniformly located on approximately parallel lines normal to the missile tra- 
jectory. The_area per ship in the BG is clearly ar ) which can be visualized as a square 
Sf side, ( a) , See Wicure 3.a, Hence the approximate distance ( D ) between ships 


1S 


= /#4. 

D= No 
As seen in Figure 3.b, the random variable L is the distance which the missile travels as 
it traverses the BG. The number ( N’) of possible ship encounters by the missile can be 


approximated as the integer value closest to N where, 
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Figure 3. Missile attacks and locations of the ships in BG for the model 


Lower case n, |, and, n’ will denote particular realizations of the r.v.’s N, L, and, N’ Since 
2R is the maximum distance which a missile can travel within the BG, the maximum 


number of encounters will occur for L = 2R; this 1s given by 


No 
Ivers =" 21K 7 


We assume the missiles meet at least one line of ships in the BG. In other words we are 


not concerned about missiles which miss the BG. Hence the range of N 1s, 


l<N<2R - ; 
A 
In the absence of intelligence reports and any detection of enemy missiles, a missile at- 
tack may come from any direction. Since there is a circular symmetry because of random 
direction of attacks, we assume the BG has a circular shaped formation. A circular for- 
mation can give more laver defense to the high-value ships in the BG by locating them 
in or near the center in case of missile attacks from any direction, and it can give efficient 
ultra high frequency ( UHF )/ visible communication conditions among all ships in the 
BG and the BG commander. 


In a circular formation, the area ( A ) covered by the BG 1s 
A=nR’. 


So the maximum number of possible encounters is 








No No 
Nmax = 2K ee cs (daly) 
mR 
and for any distance L 
Net |= (3.2) 
zR? | 


Let the closest point of approach ( CPA ) to the center of the BG for an ASCM be 
the random variable, Y, which is uniformly distributed in |0,R| From Figure 4 we see 
that Y, determines L, and therefore also N. 

Let N’ be a random variable which is N rounded to the nearest integer. Then we 
can find its probability mass function values ( PMF ) based on the underlying uniform 
distribution of Y, for any given N, { the number of ships in the BG ); see Figure 4. 





Figure 4. Drawing for calculation of the eqn. 3.2 and 3.3 


From Figure 4 we see that 
L=2./R’-Y’ , (3.3) 


hence, using eqns. (3.1), (3.2), and, (3.3) 


9) pe vee 
Naf eas ee = ie (marae 


If we take the square of the each side, 


2 4 (NoR?= NY) AN, 4N,Y° 
\ = 
R TLR 
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Solving for Y we find, 


~2 nN? 
- ( 4No ); 


N°x 
Y=R,/(1 ~ (ax) 


Replacing the r.v.’s Y and N with y and n we obtain the correspanding relations between 


4No 
Tl 


Pe 


the realizations of the r.v.’s. The latter equations are used to obtain the derivative in the 


transformation formula 


fal) =f IL , 


hence, using f/,{y) = 7 


and finally 


(3.4) 





Part of the computer program which we discuss in Chapter 4 computes the proba- 
bility mass values ( PMF ) for a discrete number of encounters using the PDF of n which 
was calculated in eqn. (3.4). As seen in Figure 5, the dotted line shows the points which 
are calculated by using eqn. ( 3.4 ) for N, = 14. Since there are a discrete number of 
lavers in reality we find the probability mass values ( PMF ) for each discrete layer by 
calculating the sub portions of area under the dotted line which are; 0.5 to 1.5 for Ist. 
layer, 1.5 to 2.5 for 2nd. laver, 2.5 to 3.5 for 3rd. layer. For Figure 5 the maximum 
number of layers = 4.2 which is calculated using the eqn. ( 3.1 ) for Nj = 14 1n the BG. 
Since the maximum N’ 1s the integer number closest to n, which is 4, we then use 3.5 to 
4.2 for the probability of the 4th. layer. 

Now let us name the small portions of the area under f,{n) as, A, A;, A;,..., A, Where 


11S the closest integer vaiue to ,,,. If we want to calculate the area under the dotted line 


( PDF of 'n ), which will be called P in this calculation, we can first calculate the A, which 
is the area under the /,() from 0 to 0.5, 


Vie [Jat E | en 


Then we can find the area of probability for PDF of n between 0.5 and 1,,,, by calculat- 


ing, ( since N > 1, we eliminate probability associated with N’ = 0.) 
P=1—Ap 
or we can find total area by calculating 
P=A,+A,+A3+...4+4; , 


where 


5 ee 
4\=| $$ dn , 
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Figure 5. Calculation of the PMF values for # encounters per missile 


Either of these two approaches should give the correct result for total probability under 
the PDF of n=>0.5. We shall use these two approaches to check the results of the com- 
puter program presented in Chapter 4. 

After finding the total area under the PDF of n> 0.5, we normalize probability mass 


values of each laver with; 


Ay 
PMF, =F * 

Ay 
PMI, = > 


1S 


A; 


PMF, = ae 
and for ith. interval, 
PME. =< 
i P ° 


Now the sum of all the probabilities should be 


PMF, + PMF, + PMF,+...+ PMF,=1 , 


so that the number of possible encounters is an integer > 1 rather than > 0. 
Figure 5 gives an example calculation of the normalized PMF values for the number 
of possible encounters when N, = 14. The dotted line shows the values of f,{(”) and the 


solid line connects the normalized PMF value points for each discrete layer. 
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B. THE MODEL TO CALCULATE THE PROBABILITIES OF SHIP HIT, 
MISSILE KILL, MISSILE DIVERSION, AND SHIP SURVIVABILITY 
I. Probability calculation for one layer | 
The probabilities of ship hit , missile kill, mussile diversion, and ship survivability 
can be calculated by constructing a sample probability tree model for one layer. As seen 
in the tree diagram in Figure 6, when a missile arrives to a line of possible encounters it 
may mect one or more ships or no ships. We now assume that the ships are located 


randomly along the line perpendicular to the missile trajectory. 
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Figure 6. Probability tree diagram for one layer of encounters 


As can be seen in Figure 7, if a is the angular search beam of the missile, then the 
probability a missile does not detect a ship on the layer ( PONS ) is 


PONS =e” , (3.5) 


Prine = V¥N/A 


@ P,[0 ships in missile beam-width] * exp{-2«) 


Figure 7, Calculation for probability if missile meets at least a ship 


which is based on the Poisson Approximation to the Binomial Distribution of the num- 
ber of ships falling within the angular search beam width of the missile. On the other 


hand the probability a missile detects at least one ship per layer ( POS ) is obviously 
POS=1-—e7** . (3.6) 


a. Missile meets at least one ship 
If a missile meets at least one ship, then it can be successfully distracted or 


not distracted by the ship using a distraction decoy with quality factor Q. The dis- 


traction probability ( POD ), which is the ratio of the equivalent number of false targets 
( QD ) to the total number of targets (1 + QD )I , and is assumed to be given by: 


OD 


POO= Te OD) 


(3.7) 
where D is the number of distraction decoys deployed by a given ship during the missile 
attack. On the other hand, the probability of no distraction of a missile given it meets a 
ship ( POND ) is 


POND = ( (3.8) 


eae eee 
ea Oe) 
Then: 

If distraction occurs, then the missile can be destroyed by a hard kill; The 
probability of hard kill given the missile is distracted 1s 


OD 


—_ — —2a ——— SD 
PMK,=(1-e ~)( (1+ OD) 


meus (3.9) 
If the missile is not hard killed it can splash down and be destroyed ( soft kill ) or not 
splash down and instead be diverted after the distraction of the missile occurs. The 
probability of splash down of a missile = PSPL. If splash occurs obviously we can ac- 
cept this as a missile soft kill. On the other hand, even if it is not splashed it is still di- 


verted. In this case the probability of missile kill given the missile is distracted is 


OD 


— _——_ —2a AD 
PMK,=(l—e ~)( (1+ OD) 


\(1— PHK)PSPL , (3.10) 


and the probability of diversion ( PDIV ) given the missile is distracted is 


- OD 
V, =(1l-—e "(= (1 - ~ 
PDIV,=(1—-e ~)( (1+ OD) \ — PHK)(1 — PSPL) (nel) 
If a missile is not distracted with the probability calculated in eqn. (3.7) it 
can be seduced or not seduced by the BG using seduction decoys. In this case the 


probability of seduction ( PSED ) is 


1 Center for Naval Analysis ( CNA ) 86-1737 / 15 September 1986 


PSED =(1—(1—B)”>) , (3.12) 
and, the probability of no seduction 1s 
PONS =(1—-— B)”” (3.13) 


respectively, where 
1. B = Probability of breaking the lock of a missile for a seduction decoy in the BG. 


2. MS = Number of seduction decoys in per ship. 


If seduction occurs the missile can be destroyed by hard kill or not destroyed by hard 


kill. The probability of missile is hard killed given the seduction occurs 1s 


l 


— _ —22 Sa 
PMK,;=(l-e ~)( (1+ OD) 


(1 — (1 — B)?)PHK . (3.14) 
The probability of no hard kill given that seduction occurs is 


PMK, =(1—e 7”)( (i =(1 = By) eee (3.15) 


a 
(1 + QD) 
and the probabability of diversion given the seduction occurs is 


PDIV, =(1—e 7°) (t= = By yer Sese7) (3.16) 


l 
(1 + QD) 
If seduction does not occur, then the probability of hard kill given the seduction does 
not occur 1s 


PMK, =(1—e *”)( B)”°)PHK , (32m 


l 
—<——$ i 
or if it cannot be hard killed it may hit a ship so the probability of a missile hitting a ship 


given the seduction does not occur 1s 


l 


(1 — eta —_1 
PSHI =(1-e"™\(7 op 


(1 — B)°\(1 — PHR) . (3.18) 
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b. Missile meets no ship 


If a missile meets no ships, then it can still be destroyed by hard kill, then 


the probability of hard kill given the missile meets no ships 1s 
PMK, =e "*PHK , (3.19) 
and the probability of diversion given the missile meets no ships is 
PDIV, = e~*"(1 — PHK) . (3.20) 


From all of these equations we calculated above, we can calculate the total 


missile kill probability ( TPMK ) per one line of encounters by using the equations. 


TPMK = PMK, + PMK, + PMK, + PMK, + PMKs + PMKz . (3.21) 


If we write the whole equation for TPMK; 


—2% OD —24 OD 
TPMK=(l—e MT Ob) PAA FU —e XT OD) — PHK)PSPL 
—22 l MS —22 l 
+(l-e “)( (1 + OD) \l —(l—- B)"~)PHK+(l—e ~)( (1+ OD)? 


1-1 — BY) — 2a l __ pS 
(l1—(1—8)°~)(1 — PHK)PSPL+(l-—e “) (1+ OD) (1 — B)"~)PHK 
+(e -")PHK . (3.22) 
The probability of ship hit in the BG ( PSHI ) is again, 


po (ee ——— > (1 — BY’) — PHR) , (3.23) 


(1+ 0D 


so survivability of BG ( the probability that all ships survive the attack ) as a measure 
of effectiveness for one line ( PSSUR ) is 


_ x. .. = _ — pyMS;y _ 
PSSUR = (I~ (1 e™")( 7 py MM — BY" — PHR)) (3.24) 


The total diversion probability ( TPDIV ) is then 
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TPDIV = PDIV, + Pei oa 


OF 


fn, —2a OD —2a l 
TPDIV =(l—e (ToD + OD) \(1 — PHK)(1 — PSPL) + (1-e (Ts oD TOD) ) 
(1 — (1 — B)””)@ — PS PE) es Pike, (3.25) 


Equations 3.21, 3.22, 3.23, 3.24 are calculated on the probability tree dia- 
gram for one line of encounters in Appendix B. The reader can refer to Appendix B to 
get an overall idea for understanding better. 


2. Probability calculation for all the layers 

We have shown how to calculate the maximum number of lines using eqn. (3.2), 
the probabilities of any ship being hit in the BG, missile kill, missile divert, and 
survivability, for one line of encounters. Now we can find the probabilities conditioned 
on the number of layers ( or potential encounters with ships ) which make up the whole 
BG. For these calculations we shall use a probability tree, which can be seen in Figure 
8, for all the lavers in the BG, and we will be using the equations that we calculated for 
one layer encounters. 

As we can see in Figure 7, for each line the missile may hit a ship or get diverted 
or get ” killed ”. Ifa missile hits a ship or gets killed, the threat of the missile is termi- 
nated. But if a diversion occurs the same conditions and probability calculations take 
place again and again until it passes the last line of encounters within the BG. For the 
last line, diversion is equivalent to the termination of the missile. 

By using the equations in Part | and also using Figure 7 and Figure 8 we can 
find the following probabilities for a BG: 

Probability of missile destruction for a BG ({ PMKBG ) averaged over the 
number of lines of encounters for a given number { N ) of ships is 


n 
PMKBG = ) Pr[MK | i layers} PMF, . (3.26) 


i=] 


For instance by adding the probabilities for mutually exclusive ( M.E.) events we have: 
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Figure 8. ‘Probability tree diagram for all the layers in the BG 


Pi{MK |3 layers]= PMK + (PDIV) PMK + (PDIV’) PMK . 


Probability of a ship getting hit for a BG ( PSHBG ) which has n lines of en- 


counters for a given number ( N ) of ships ts 


nN 
PSHBG = ) Pr [SHI |i layers] PMF, . (3.27) 
i=] 


For instance by adding probabilities for M.E. events we have, 


Pr[SHI |3 layers] = PSHI + (PDIV) PSH + (PDIV’) PSH . 


2\ 


Probability of diversion of a missile for a BG ( PDIVBG ) which has n lines of 


encounters for given number ( N ) of ships is 


PDIVBG = ) (PDIV) PMF, . (3.28) 
=) 


Also we can calculate the probability of ship survivability ( PSSBG ) as measure 
of effectiveness ( MOE ) for the BG. 


PSSBG =(1 — PSHBG) . (3.29) 


We will be using the probabilities we calculate in eqns. 3.26, 3.27, 3.28 and 3.29 
for drawing the graphs to understand the different relationships between the survivability 
of the BG and factors such as the number of each different type of decoy, the number 
of ships in the BG, the number of missiles attacking the BG, the distraction decoy Q 
factor, the seduction decoy B, factor, etc. 

To illustrate these relationships we can calculate probabilities in eqns .3.26, 3.27, 
3.28 and 3.29 by using the FORTRAN 77 program in the Appendix. After we run the 
program by changing the different parameters we get the example results as Figures 9 - 
17. 

As seen in Figure 9 we run the program with the parameters: 
Number of ships in the BG (N,) = 5...50. 
Number of distraction decoys per ship ( DI ) = 1. 
Number of seduction decoys per ship ( MS ) = 1. 
Hard kill probability of BG ( PHK ) = 0.25. 
B factor of seduction decoys = 0.2. 
Q factor of distraction decoys = 0.2. 
Splash probability of missile ( PSPL ) = 0.05. 
Alpha angle of missile (a) = 15”. 


ae ee ee ee lS 


and later we get the ship survival probabilities as the number of ships in the BG are in- 
creasing. We can easily see from the graph that when the number of ships increases in 
the BG, survivability of ships decreases. Also ship survivability of BG has a more nega- 
tive slope ( decreasing survivability ) for N, = 5 to 15 than N, = 15 to 50 ( 50 can be 
thought of as very large = infinite number of ships for BG ). When the number of ships 


goes to infinity, then the probability of survivability of BG is almost constant. 
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Figure 9. =‘ Probability of no hits, while number of ships varies 


If we look at Figure 10 and the Figure 11 we can easily see that the different 

types of decoys affect the survival probability differently. In Figure 10: 

1. The number of ships (N,) varies from 5-30. 

2. Number of distraction decoys ( DI ) varies from 1-5. 
We get five different survival curves for BG. while the other parameters remained un- 
changed and are, 

3. Number of seduction decoys per ship ( MS ) = I. 

4. Probability of hard kill level of BG ( PHK ) = 0.25 

5. Splash probability of the missile ( PSPL ) = 0.05 


6. B factor of seduction decoys = 0.2 
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SURVIVABILITY OF BG 
N=5S-30, Di=1—-5, MS=1, PHK=.25, PSPL=.05, B=Q=.2 
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Figure 10. =‘ Probability of no hits, distraction decoys vary from | to 5 


7. Q factor of distraction decoys = 0.2 

8. Alpha angle of the missile (a) = 15° 
In Figure 11: 

1. Number of distraction decoys ( DI ) = 1. 

2. The number of seduction decoys ( MS ) varies from 1-5. 
We get five different survival curves for the BG. while the other parameters remained 
unchanged and are the same as above. As we see in Figures 10 and 11, different types 
of decovs affect the survivability. In this case the number of seduction decoys has 


slightly more effect than the distraction decoys on survivability of a BG. 
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Figure 11. Probability of no hits, seduction decoys vary from 1 to 5 


If we change the number of decoys, but we do not change the number of ships 
in the BG, we can get Figure 12 and Figure 13. In Figure 12: 

1. Number of distraction decoys ( DI ) varies from 1-5. 

2. The number of ships (N,) = 5. 

3. Number of seduction decoys per ship ( MS ) = 1. 

4. Probability of hard kill level of BG ( PHK ) = 0.25. 

5. Splash probability of the missile ( PSPL ) = 0.05. 

6. B factor of seduction decoys = 0.2. 

7. Q factor of distraction decoys = 0.2. 

8 


. Alpha angle of the missile (a) = 15’. 
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Figure 12. Survivability of BG while num.of distraction decoys varies 


In Figure 13, we modifv the first and third parameters so that: 
1. Number of distraction decoys ( DI ) = I. 
3. Number of seduction decoys ( MS ) varies from 1 to 5. 
In these two figures, when the seduction decoys were changed, we obtain a - slightly 
higher survival probability compared to when the number of distraction decoys were 
changed. The benefit of expendables are apparent from these figures. 
The number of missiles attacking the BG is also an important factor for 


survivability, which is: 


Pr ( Survival given M missile attacks ) = Pr{ Survivability of BG ye , M=1,2.... 
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Figure 13. Survivability of BG while num.of seduction decoys varies 


This assumes independence for each missile and that the decoys are not depleted during 


successive attacks. After running the program and calculating the survival probabilities 


for M missiles we got the graphs in Figure 14 and Figure 15. In Figure 14 we run the 


program for M missile attacks while, 


l. 
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Number of missile attacking the BG ( M ) varies from | to 10. 


. Number of distraction decoys ( DI ) varies 1 to 5. 

. Number of seduction decoys ( MS ) = 1. 

. The number of ships (N,) = 10. 

. Probability of hard kill level of BG ( PHK ) = 0.70. 
. Splash probability of the missile ( PSPL ) = 0.05. 
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Figure 14. Survivability of BG while num.of distraction decoys varies 


7. B factor of seduction decoys = 0.5. 
8. Q factor of distraction decoys = 0.7. 
9. Alpha angle of the missile («) = 15°. 


In Figure 15 we run the program for M missile attacks making the following changes: 
1. Number of seduction decoys ( DI ) varies | to 5. 


2. Number of distraction decoys ( MS ) = 1. 


From Figure 15 and 16, we can see that seduction decoys have more effect on 
survivability than distraction decoys have. Especially after the third decoy has been de- 
ployed. For each type of decoy there is not a very big difference on survivability for 
number of attacking missilies ( M ) = 1, but when M increases the seduction decoys 
have more effect yielding significantly higher survivability. 

The specific values of Q and B are important factors for determining the effect 


of expendibles. Figures 16 and 17 show the effects of varying the Q factor of a dis- 
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Figure 15. Survivability of BG while num.of seduction decoys varies 


traction decoy and the B factor of a seduction decoy, for the case: | seduction and | 
distraction expended for engagement by | missile. In each case B and Q varied between 
0.0 to 1.0 while the other parameters remained unchanged. The effect of seduction decoy 
B factor to survivability is higher than the effect of distraction decov Q factor. In surface 
ships Q and B may assume any value between 0 to I. The specific values of B and Q 
depend on systems employed and the dynamics of the engagement. In Figure 16, we run 
the program for different values of distraction decoy Q factor while other parameters 
remain unchanged. 

Q factor of distraction decoys varies from 0 to 1. 

Number of seduction decoys ( MS ) = I. 

The number of ships (N,) = 15. 

Number of distraction decoys ( DI ) = 1. 

Probability of hard kill level of BG ( PHK ) = 0.25. 

Splash probability of the missile ( PSPL ) = 0.0S. 


DnB wn o 
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Figure 16. Survivability of BG while Q factor of distraction decoys varies 


7. B factor of seduction decoys = 0.2. 


8. Alpha angle of the missile (a) = 15°. 


In Figure 17 we run the program for different values of the seduction decoy B factor 


while other parameters stay unchanged. 
1. B factor of seduction decoys varies from 0 to 1 ( 0.1 increments ). 


2. Q factor of distraction decoys = 0.2. 
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Figure 17. Survivability of BG while B factor of seduction decoys varies 
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IV. CONCLUSION 


While defending the BG against anti-ship cruise missile attacks in the open ocean, 
a major concern is survivability of the BG and the impact of the different performance 
factors on the survivability of the BG. 

Since there 1s a big difference between the speed factors of a BG and a missile, we 
have assumed that the BG’s location is fixed in this particular situation. We also defined 
the measure of effectiveness ( MOE ) as the probability that all ships survive the attack 
of M missiles. 

Major factors whiel impact the BG’s survivability are number of ships in the BG, 
the BG’s hard kill level, and the BG’s expendable decoys. We have seen the benefit of 
expendables in the examples of the previous chapter which were drawn using the data 
calculated by the FORTRAN 77 computer program in the Appendix. Specifically, we 
have investigated the sensitivity to the distraction decoys Q factor and the seduction 
decoys B factor which are important for determining the effects of expendables. 

The analvtic results of the model which were calculated by the computer program 
give results consistent with those expected in the beginning. If the number of ships in- 
creases and the other factors remain constant, survivability of BG / ship decreases, in 
part because of the existence of the nonlinear interactions between ships in the BG. The 
tactics of single ship defense against ASCM’s, the position of the ship in the BG ( the 
ship’s nearness to center / high-value ships ) may also affect the BG’s multiple ship 
overall survivability, or the survivability of high value targets such as CV's. 

Different types of decoys affect survivability differently. This 1s done by either the 
B or Q values of the different type of decoys, and / or the number of decoys that can 


be deployed. Specifically distraction decoys’ Q factor and seduction decoys’ B factor are 
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important for determining the effects of expendables. Hence one should analyze 
survivability of the BG and include cost factors of the expendables in the model before 
allocating and deploying them on each ship. Different trade-offs also can be made be- 
tween decoys, hard kill levels, and number / position of ships in the BG, which has lim- 
ited space and budget for devices, weapons, amount of decoys, and missiles. 

When we include the cost factors of these variables in our model we can easily get 
trade-offs by using a survivability cost comparision between different types of assets. 
By using the results of this model we can change the different variables such as the 
number of decoys and hard kill levels, subject to given cost constraints, to optimize the 


survivability. 
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APPENDIX A. FORTRAN PROGRAM TO CALCULATE THE SHIP 
SURVIVABILITY IN BG 


TeFedeFe Hse HR RRR KEK RRIF REEREEKEEKERRRERRRRRE EKER RERERERERRE KK 


* 21,AUG,1988 * 

* PROGRAM SENGEL MOD7 * 
* AUTHOR. CENGIZ SENGEL * 
* NPGS OPERATIONS RESEARCH. ( 360 ) * 
% Kk KK ¥ 
¥ THIS FORTRAN 77 PROGRAM CALCULATES THE SHIP SURVIVAL PROBABILITY * 
* IN THE BATTLE GROUP ( BG ) FOR N GIVEN NUMBER OF SHIPS AND OTHER * 
* PARAMETERS WRITTEN BELOW. PROGRAM ALSO CALCULATES THE SHIP * 
* HIT FOR BG, SINGLE SHIP HIT / SURVIVAL, ATTACKING MISSILE, * 
* DIVERSION, KILL PROBABILITIES FOR M MISSILE ATTACKS. ASSUMING * 
* INDEPENDENT MISSILE ATTACKS M CAN BE 1,2,...,N AND BG HAS ALMOST * 
* CIRCULAR FORMATION AGAINST MISSILE THREAT WHICH MAY COME FROM * 
* ANY DIRECTION. ve 
* KIekvedededewde Foxe ye AY 
* PARAMETERS: * 
* N= # OF SHIPS/ INPUT * 
* L=MAXIMUM LAYERS FOR GIVEN N CREAL) * 
* MXL= L VALUE ROUND TO NEAREST INTEGER % 
* TSUM= TOTAL PROBABILITIES BETWEEN .5 AND L * 
* TS = TOTAL PROBABILITIES WITH USING 1-PROBABILITY BETWEEN 0 * 
* AND .5 * 
* TSUM= TOTAL PROBABILITY AREA UNDER PDF FCN BETWEEN 0.5 ,MXL * 
% A=ALPHA ANGLE OF ONCOMING MISSILE/ INPUT * 
* Q= QUALITY FACTOR OF DISTRACTION DECOY (Q. LE. 1)/INPUT % 
* D= # OF DISTRACTION DECOYS/ INPUT * 
* B= PROBABILITY OF BREAKING THE LOCK FOR A SEDUC. DECOY/INPUT * 
* MS= # OF SEDUCTION DECOYS/ INPUT * 
¥ PHK= HARD KILL PROBABILITY OF SHIP FOR KILLING MISSILE/INPUT* 
7 PSPL=PROBABILITY OF MISSILE GET SPLASHED AFTER DECOYED/INPUT * 
* PSH= PROBABILITY OF SHIP GET HIT * 
% PDV= PROBABILITY OF MISSILE GET DIVERTED * 
% PMK= PROBABILITY OF MISSILE GET KILLED ( HARD KILL OR SPLASH)* 
% CPSS= CONDITIONAL PROBABILITY OF SURVIVAL FOR BG ( PROB. OF * 
** NO HIT) * 
* CPDV= CONDITIONAL PROBABILITY OF MISSILE DIVERSION * 
* CPMK= CONDITIONAL PROBABILITY OF MISSILE KILL * 
7 CPSS1= CONDITIONAL PROBABILITY OF SINGLE SHIP SURVIVABILITY * 
IN THE BG * 
CPSH= CONDITIONAL PROBABILITY OF SHIP HIT IN THE BG * 


CPSH1= CONDITIONAL PROBABILITY OF SINGLE SHIP HIT alee = 
IN THE BG 
CPSSM= CONDITIONAL PROBABILITY OF SHIP SURVIVABILITY IN Cao * 


OF M MISSILE ATTACKS 
devededeeedeved kere eK * 


Setetetestededetede fete sedetede te fe deteTehe TREAT RRRRRIRRERERERRERRERRITRRERERRERERERREERERERRRE 
REAL N,L,MXL,P,PRO1,PL,MXL1,P0,P1I,TSUM,TS,A,Q,D,B,MS, PSPL, PHK, PSHI 
5 ,PSH, PDV,M,PDVI, PMKI, PMK, DI ,MXL2,PL1, PLO,AL, CPSS,CPSH,CPSS1,CPSH1, 


+b tb ye tt 
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6CPSS1M,CPSH1M,K,R 

CALL EXCMS ('FILEDEF 01 DISK DEMET OUTPUT Al1') 

DATA N/10/,A/15. /,PSPL/. 05/,PHK/. 7/,DI/1. /,Q/.7/,B/. 5/ 
* DEFINING THE CONSTANT PI AND ALPHA RADIAN 

PI=3. 141592654 

AL=-2*(A*(PI/180. )) 


DO 14 R=1. ,5. 
* DO 14 N=5., 50. 
M=N 
MS=R 
* WRITE(01,24) N 


24 FORMAT(1X,F5. 1) 
*CALCULATING THE MAXIMUM NUMBER OF LAYERS 
L=2*SQRT(N/PI) 
*ROUNDING THE MAXIMUM NUMBER OF LAYERS TO THE NEAREST INTEGER 
MXL=ANINT(L) 
SUM=0 
MXL1=(MXL-1)+. 5 
*CALCULATION OF THE AREA UNDER PDF=TS (.5,L) AND AREA OF LAST PORTION 
PO=ABS(1-SQRT(1-(. 5**2*PI/(4*N)))) 
PL=(SQRT(1. -(L*2. *PI/(4. *N))))-(SQRT(1. -(MXL1**2. *PI/(4. *N)))) 
TS=1-P0 
PLP=ABS(PL)/TS 
cs WRITE(01,22) PLP,L,TS 
22  FORMAT( 1X,F9.8,1X,F9.6,1X,F9. 8) 
T=0,0 
CPSH=0. 0 
CPMK=0. 0 
CPDV=0. 0 
PSH=0. 0 
PMKI=0. 0 
PDV=0. 0 
*CALCULATING THE PROBABILITIES FOR ONE LAYER YOU CAN FIND PMK BY TWO WAY 
PSHI=(1-EXP(AL) )*((1/(1+(Q*DI) ) )*(( 1-B)**MS )**( 1-PHK) ) 
PDIV=( 1-EXP( AL) )*(((Q*DI) /(1+(Q*DI)) )*( 1-PHK)*( 1-PSPL) 
1+((1-(1-B)***MS)*( 1-PHK)*(1-PSPL)*(1/(1+(Q*DI)))))+ 
2( (EXP( AL) )**( 1-PHK)) 
*FIRST WAY OF FINDING PMK 
* PMK=1-(PDIV+PSHI) 
*OR SECOND (LONGER WAY) OF FINDING PMK 
PMK1=( (1-EXP(AL))*( (Q*DI) /(1+(Q*DI) ) )*( PHK) ) 
PMK2=( ( 1-EXP( AL) )*((Q*DI) /( 1+(Q*DI) ) )*( 1-PHK)*PSPL) 
PMK3=( 1-EXP(AL) )*(1/( 1+(Q*DI)))*(1-(1-B)**MS )*PHK 
PMK4=(( 1-EXP( AL) )*(1/(1+(Q*DI) ) )*( 1-( 1-B)**MS)*(1-PHK)*PSPL) 
PMK5=( ( 1-EXP( AL) )*(1/(1+(Q*DI) ) )*((1-B)**MS )*PHK)+EXP(AL)*PHK 
PMK6=EXP( AL) *PHK 
PMK=PMK 1+PMK2+PMK3+PMK4+PMKS5 
*CALCULATION OF PROBABILITIES FOR LAYERS 
DO 15 F=1,MXL,1 
PR=0. 0 
PSH=(PSH+PDIV**(T)*PSHI) 
T=T+1 
PDV=PDIV**F 
IF(F. EQ. MXL)THEN 
CPSH=CPSH+PSH*PLP 
CPMK=CPMK+PMKI*PLP 
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CPDV=CPDV+PDV*PLP 

ELSE 

PR=SQRT( 1-((F+. 5)%**2*PI/(4*N) ))-SQRT(1-((F-. 5)**2*PI/(4%*N))) 
PRO1=ABS(PR)/TS 

PMKI=(1-(PDV+PSH) ) 

CPSH=CPSH+PSH*PRO1 

CPMK=CPMK+PMKI*PRO1 

CPDV=CPDV+PDV*PRO1 


END IF 
WRITE(01,27)DI,MS 
7 FORMAT( 1X,2(F9. 6,1X)) 
15 CONTINUE 
CPSS=1-CPSH 
*DISPLAY CONDITIONAL PROBABILITY OF SHIP HIT, MISSILE KILL, DIVERSION 
* WRITE(01,23)CPSH,CPMK,CPDV 


*DISPLAY THE SHIP SURVIVABILITY AS 'MOE' 
*SINGLE SHIP PROBABILITY CALCULATIONS 


CPMK1=CPMK/N 
* CPSH1=CPSH/N 
* CPSS1=CPSS/N 


*BG PROBABILITY CALCULATIONS GIVEN M MISSILE ATTACKS TO THE BG 
CPSHM=1. -( 1. -(CPSH)**M) 
CPSSM=1. -(1. -(CPSS)**M) 

*SINGLE SHIP PROBABILITY CALCULATIONS GIVEN M MISSILE ATTACKS TO THE BG 

*WHICH IS M= 1, 2, 3, ...,; N. 
CPSH1M=1. -(1.-(CPSH/M)**M) —- 
CPSS1M=1. -(1. -(CPSS/M)*#*M) 
WRITE(01,23)N,MS,DI,CPSS,CPSSM 

23.  FORMAT( 1X,F4.1,1X,F4.1,1X,F4. 1,1X,2(F9. 7,1X)) 

14 CONTINUE 

31 STOP 
END 
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APPENDIX B. PROBABILITY CALCULATIONS FOR ONE LINE OF 
ENCOUNTERS. 


We can obtain all the equations calculated for one line of encounters with the aid 
of the probabilities on the probability tree diagram in Figure 18. 


For calculation of the total missile kill probability ( TPMK ): 


= —2a OD 
PMK, =(1—e "(7 oy PAK , 


| = —2a Ob 
PMK, =(1-e°2*\_———— T+ OD) \(1 — PHK)PSPL , 


l 


S le, es 2 MS 
(1+ OD) (l= 3) SPH , 


PMK, =(1 —e7*"\( 


PMK,=(1-e 1-—(1—B)”°)\(1 — PHK)PSPL , 


I 


ees a MS 
Te OD) Mi - BY PHK 


PMK; =(1—e 7”)( 


PMK, =e “PHK , 


TPMK = PMK, + PMK, + PMK; + PMK, + PMKs + PMK, 


_ ale OD 7 eu LaeD _ 
TPMK =(l—e (Toby OD) )PHK+(1-e 2%) ——=——. T+ OD) \(1 — PHK)PSPL 


ah resre -)( d= C= Beek Pte 7" 


] 
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oF 
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eel B) >) =PHR)PSPL + (1—e% “\ yi — BY’? PHK 


(1+ QD) 
+(e°")PHK . 


For calculation of the total missile diversion probability ( TPDIV ): 


« -2,_ QD, x 
PDIV, = (1— e)( 7 oa \( — PHK\( — PSPL) , 


PDIV, = (1 — e 2%) ial Be )\dueS PL). « 


1 
cals) 
PDIV, =e °*(1 — PHK) , 
TPDIV = PDIV, + PDIV, + PDIV, . 


l 


aa —2% OD mae Te oe 
TPDIV =(1 — 7 \1 — PHK)(1 — PSPL) + (1 -—e"*”)( (1 + OD) 


(I+ OD) Bs 
(ed=5))\(1 — PSPL) +(e) — PH) . 


For calculation of the ship hit probability and ship survival probability: 


l 


it gate P= PH®) 


PSHI=(1-—e° 7") 


Bee ey SS 
PSSUR=(1-(l-e “)( (1+ OD) )\1 — By)" (1 — PHK)) . 
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